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ABSTRACT 



v„ 

In this thesis, an interactive computer program for the 
analysis and design of time invariant unity feedback control 
systems is presented, using cascade or feedback or both 
types of compensation. 

By using this program, the user is freed from the 
tedious, time consuming and error prone method of hand 
calculations, letting the computer handle these tasks 
efficiently and speedily. The user can then concentrate 
fully on the placement of poles and zeroes of the 
compensator ( s) used. 

Design of control systems by classical methods being 
essentially a repetitive, trial and error procedure, this 
program greatly cuts down the turn around time and leads to 
faster, more satisfactory results. 
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INTRODUCTION 



I . 

During the past two decades, the scientific and 
engineering coiximunities have witnessed an ever increasing 
role of the digital computer in the fields of research, 
development and analysis of systems. The computer, today, 
is being used to solve engineering problems, whose solution, 
until very recently required long and tedious procedures. 
Still, it is probably true that this machine has potential 
not fully recognized yet, which is why the attention of so 
many computer scientists /engineers is focused on devising 
more efficient and innovative operating procedures. 

Control system design is one area where classical 
theory has been extensively developed and used. It is fair 
to say that even today, most analysis and design problems of 
linear, time invariant control systems can still be 
approached using the methods developed by Bode, Nyquist and 
others . 

A totally new approach to the design of control systems 
became available with the development of optimal control 
theory and the state variable analysis. These methods have 
been intensively developed in the last 10-15 years, but now 
their weaknesses have been exposed too. The "states" of the 
plant may not necessarily represent physically measurable 
quantities, and consequently it may not be possible to 
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implement the results at all. Luenberger's observers, 
designed to overcome this problem, can at best provide 
estimates of the state trajectory. Furthermore, the optimal 
control approach to design relies very heavily on mathe- 
matical manipulation, providing little insight to. the actual 
working of the plant; the only -input of the designer being 
the form of the cost function. 

An intelligent use of the speed and information 
processing ability of the digital computer, coupled with the 
reliable features of classical theory appear to be the best 
solution to the problem at hand. The classical approach to 
design, being essentially a trial and error method, if the 
order of the system is fairly high, the number of repetitive 
calculations and the time required to perform these calcu- 
lations becomes prohibitively large, the assistance of the 
computer in such problems becomes indispensable. 

The work in this thesis was to develop an interactive, 
user oriented computer program that would prompt the user 
to input the transfer function and cascade/f eedback compen- 
sators. The program would then display on the IBM 3277 - 
Tektronix 618 dial screen terminals the Bode Plot of both 
magnitude and phase. The program could be repeatedly used, 
with the user having the option to change/modify the compen- 
sators, each time viewing the effect of his modifications on 
the screen until he arrives at a satisfactory solution. 
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II. CONTROL ENGINEERING ANALYSIS 



A. GENERAL 

A continuous time control system may be represented in 
one of the following forms: 

a. Transfer functions 

b. State equations 

c. System block diagrams or signal flow graphs 
Algorithms exist in almost every undergraduate control 

engineering text to convert the system representation from 
one given form to another. Gianniotis (Ref. 1) describes a 
simple method of converting from transfer function to state 
variable form in Chapter II. 

The transfer function representation in its most general 
form is: 



A S^ + A ,S™ ^ + A .s"' ^ + A S‘ 

m m-'l m-2 o 



B s”^ + B ^ + B „S^ ^ + B S*^ 

n n-1 n-2 o 



Usually the mathematical description of the system is 
found in the transfer function form in the literature. 
Analysis and design of control systems by classical methods 
also requires the representation of the system in this form. 
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This thesis does not address the problem of converting 
from one form of representation to another. It is assumed 
that any conversions necessary have already been performed 
and that the system is represented by its open loop transfer 
function . 

B . PROBLEM FORMULATION 

Each system design has its own unique characteristics, 
but in general the system has to meet some kind of perfor- 
mance standards. These performance standards are generally 
provided as numerical specifications. The first step in the 
design of a control system is to analyze the system by itself 
in the usual feedback loop configuration. This is usually 
referred to as the uncompensated system. 

Analysis of the uncompensated system almost always shows 
that the system cannot meet some or all of the given perfor- 
mance standards. Usually, additional components have to be 
inserted in the system for the purpose of altering the 
performance of the system. These components are called 
compensators. Compensation is a two step procedure, in which 
additional components (compensators) are inserted to change 
the structure of the system, and these components are then 
adjusted until the performance characteristics are satisfied. 

The theory of cascade and feedback compensation is 
discussed in detail by Thaler (Ref. 2) in Chapters 5 and 6. 
Only a brief discussion of the types of compensators is 
presented here. 
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Compensators used in the design of feedback control 
^^0 g0nerall^ classified as shown in the block 
diagrams of figure 1- 



COMPENSATOR 



PLANT 




O) CASCADE COMPENSATION 




b) FEEDBACK COMPENSATION 

PLANT plant 




C) COMBINED cascade AND F E EDB ACK C OMP E NS AT I 0 N 



Figure 1. Classification of Compensation Structures 
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Cascade compensation may further be classified into two 
further types; a high pass filter usually called a phase 
lead compensator, and a low pass or phase lag compensator. 

The selection of the type of compensator (s) to be used 
depends on a number of factors, the important ones being 
experience, personal preferences, availability, system cons- 
traints, etc. Unfortunately, there are no mathematical 
techniques to help in this selection process. Generally, 
one has to complete several designs and then choose the most 
appropriate one. 

C. COMPUTER AIDED DESIGN 

Of all the classical design methods available, the Bode 
diagram technique is generally considered to be the simplest, 
at the same time providing the most insight into system 
performance and behavior. The Bode design method may be used 
equally effectively both with cascade and feedback compen- 
sation schemes. 

The computer program developed in this thesis displays 
on the terminal, initially, the Bode diagram (magnitude and 
phase) of the uncompensated open loop system. The user may 
then select the type of compensator to be used, and feed this 
information to the computer. The display changes, now showing 
the compensated system Bode Plot. This procedure may be 
repeated iteratively, with the computer updating the Bode 
Plot with each change made. 
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It was found on the average that turn around time for a 
typical third order system with two cascade compensators is 
less than five minutes. 
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III. PROGRAMMING CONSIDERATIONS 



A. GENERAL 

Any interactive software package, such as the one 
developed in this thesis must provide a simple yet unambigous 
means of data input and output. The data must be easy to 
interpret and apply to the problem at hand. Programs 
producing highly satisfactory results but requiring long 
studying time and/or special programming skills are of 
limited use only. 

The intent of this thesis is to present such a program, 
Computer Aided Design of Linear Systems. Special care has 
been taken to develop the program so that the user has to 
invest very little time learning to use it. 

B. MAIN FEATURES OF THE PROGRATl 

The development of a user oriented interactive computer 
program in solving engineering problems requires a consider- 
able amount of programming work, contributing significantly 
to its complexity and size. 

The computer program, hereafter referred to as BODPLT 
consists of a main program, a number of programmer-composed 
subroutines, a few library functions/subroutines, and various 
subroutines of the DISSPLA graphics package. The entire 
program is written in the FORTRAN IV language. 
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In brief, the whole BODPLT package works as follows: 

A user, logged into the VM/CMS environment of the system from 
the dual screen terminals, issues the command DISSPLA BODPLT. 
The package then assumes control. The program begins its 
execution by interrogating the user and calling the appro- 
priate subroutine accordingly. All programmer-composed 
subprograms are included in the main program titled BODPLT 
FORTRAN . 

The BODPLT program has the following important features: 

- runs of the VM/CMS time sharing system. 

interrogates the user in entering all problem 

specifications from the terminal. 

- can handle up to a ninth order plant transfer 
function, six first order and one second order cascade 
f i.lters . 

prompts the user to input the parameters of 
velocity feedback, acceleration feedback or approximate 
acceleration feedback as required. 

- provides the solution in tabular form on the IBM 
3277 screen and the Bode Plot on the TEK 618 terminal. 

can provide hard copy version of the problem 
specifications and tabular output by using the RECORD ON/ 
RECORD OFF execs, and of the BODE PLOT on the Tektronix 
printer where installed. 

allows problem specifications to be changed between 

runs . 
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C. PROGRAM DESCRIPTION 

The main program is the coordination center which 
controls the calling of the supporting subroutines, in order 
to input the transfer function, cascade/feedback compensators 
and other necessary information. The tabular results and the 
Bode Plot are then displayed on the two screens respectively. 

The main program as well as the accompanying subroutines 
can be found in Appendix A. They contain a sufficient number 
of comment cards to be self explanatory. 

The following is a brief description of the performance 
and purpose of the various subroutines. 

NUMER inputs the numerator of the plant transfer 
function. 

DENOM inputs the denominator of the plant transfer 
function . 

CASCAD inputs up to 6 first order lead/lag filters. 

SEC AS inputs the numerator of the second order band 

pass/band stop filter. 

SECASD inputs the denominator of the second order 
filter . 

FETCH determines the value of the radial frequency, 
w, at the origin of the x-axis. 

DECADE determines the number of decades of frequency 
to be spanned. 

FEEDER inputs the various parameters of feedback 
compensators . 
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TITLES inputs the two lines of text as headings for 
the Bode Plot. 

In addition, the main program handles the tasks of 
displaying the tabular output and the Bode Plot of the open 
and closed loop response of the system as required. 
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IV. PROGRAM PERFORMANCE INVESTIGATION 



A. GENERAL 

The program was tested by solving several linear control 
system design problems. Very satisfactory results were 
obtained in all cases with remarkable efficiency and speed. 
The only necessary condition is proper problem formulation. 
This is true for any interactive computer program. Once the 
program is used a couple of times, the user gets the 
necessary familiarity with its working. 

The example problems presented below are used to 
demonstrate the performance and capabilities of the program. 
The examples can also help the user in formulating his own 
particular problem. The examples were selected from Thaler's 
'Design of Feedback Systems' (Ref. 2) 

B . EXAMPLE PROBLEMS 

1 . Example 4.1 : A Phase Lead Network 
a. Problem Statement 

A positioning system is single loop with unity 
feedback and forward transfer function 

5.0 

G(s) = 

s (0 . 7s + 1) (0 . 3s + 1) 
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The velocity constant is not to be decreased. Design a, 
phase lead compensator which will provide less than 1,5. 

b. Solution 

The first step is to get the transfer function 
in the required form: 

5.0 

G(s) = (4.2) 

0.21S^ + l.OS^ + l.OS 

The Bode Plot of the uncompensated system drawn 
using the program is shown on figure 2 and the tabular output 
on table 1. The uncompensated system has a phase margin of 
-20 degrees. To achieve a of less than 1.5, a phase 

margin of 44 degrees or more is required. Approximately 64 
degrees of positive phase shift are therefore needed. Two 
sections of lead filter are introduced as given below: 

(S/3.0 + 1) (S/1.5 + 1) 

G = 

2 

(S/10.0 + 1) 

The compensated Bode Plot is given on figure 
3 showing a phase margin of 35 degrees. 

It may be pointed out that the final values 
for the lead filter poles and zeroes were arrived at after 
3 iterations and this design problem was solved in less than 
15 minutes , 
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Figure 2. Uncompensated Bode Plot of Example 
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Table 1. Tabular Output of Example 4.1 



FREi 



O.xOOOOQE-Oi. 
0 .ix0776E-0i 
0. 12z7i2E-0i 
0 .135936E-01 
0.l5058^E-0i 
0. i6o8iQE-0i 
0.id4785E-0i 
0.20't697E-Ui 
0 .226734E-0i 
0 .251 i88 E-Oi 
0.27o23t)E-Ci 
0 .^08240E-0i 
0 . 34x 455 E~ Ci 
0.378249E-OJ. 
U .41900a E-Oi 
0.464158E-C1 
0 .5l4175E-0i. 
0 .569580E-01 
0.630956E-C1 
0 .098 946E-0J. 
0. 774z:63E-0i 
0.857695E-01 
0 .950117 E-01 
0 .105250E + 00 
0 . ii.0 591 c +00 
0 .129155E+00 
0.143 0 72E + 00 
0 .l5a4a9c + 00 
0.17a 567 E + 00 
0.19448OE+C0 
0 .215443 5+uO 
0.236658E+00 
0 .2o4 3 7o E + 00 
0.292O64E+0J 
0.a24421E+C0 
0 .359381E + 00 
0.3981U6E+00 
0.44i004E+00 
0 .488526 E+00 
0.541 166E+00 
0 .59948 jE+ 00 
0 «o64 0al E+00 
0. 735o40£+00 
0.814911E+00 
0.602 723E+00 
0.999999E+00 
0.110775E+01 
0 . 122712E + C1 
0 .lo5 9^5E+Ux 
0. 1505a3E + 0i. 
0.1 o6 8 09E + Oi. 
0.1o47a5E+01 
0.204696E + 0i. 
0.z26754E + 0i. 
0 .251 loaE + Oi 
0.278255E+C1 
0 .a08 239 E+Ol 
0.j41454E + 0i 
0.378247E+01 
0 .419007 £+01 



MAGNITUOE 



PHASE 



0. 53979o£ + 02 
0. 530907E+02 
0. 522019E + 02 
0. 513130E+0Z 
0. 504c42E+02 
0. 495354E + 0^ 
0 . 4 864 66 E + 02 
0. 477579E+02 
0. 468o91£+0^ 
0. 459O04E+02 
0. 450918E + 02 
0. 4420a2E+02 
0. 433146E+02 
0. 4242 62E + 02 
0. 415j7oE+02 
0. 40o49oE +02 
0. 397oi5£ + 02 
0. 388 736E + 06 
0.379859E + 02 
0. J70985E+02 
0.3o2114E+02 
0.353247E + 0^ 
0. 344385E + 0^: 
0. 335529E+02 
0. 32oo80£ + 0^ 
0. 317841E + 02 
0. 309012E + 02 
0.300196E+02 
0.291397E+02 
0. 686617E + 0<^ 
0. 2738O0E+02 
0. 265130E+02 
0. <:5o43aE + 02 
0.247775E+0^ 
0. 2391O2E + 02 
0.230600E + O2 
0. 22^096c +02 
0. 2 1.60 5 5E + 

0. 205278E+02 
0. 19o96oE + 02 
0. i8a709E + 02 
0. 180467E+02 
0. 172265E + 02 
0. 163984E + 02 
0. 1555 01c +02 
0. 14o88iz + 0 ^ 
0. 137a07E+02 
0. 128181E + 02 
0. 117855E+02 
0.1 Oo 723 E +02 
0. 947399E + 01 
U.819452E + 01 
0. 6OH498E+01 
0. 544ii2E+0i 
0. 39999aE+01 
0. 253o84E+01 
0. 106322E+01 
-0.41a7o2E+00 
- 0. 1891 7o£ + Oi 
-0. 3O7Z25E+01 



-0. 90j.la4E+0i 
-C.90iz58E+0i 
-0.901a9aE+02 
-0.601546E+0Z 
-0.9U1714E+02 
-0. 90i900£+0i 
-0.902106E+02 
-0.902335E+02 
-0.902588E + 01 
-C.9028O8E+02 
-0.9J3179E+02 
-0.903524E + 02 
-0.903905E+02 
-0.9043z9E+0i 
-0.90479oE+0i 
-0.905al8E+02 
~C. 905o9oE + Ol 
-0.9 06536E + 01 
-0.907247E+01 

— 0.908036E+u2 
-0.9086 laE+02 
—O. 909887 E+ 02 
-C.610672E+02 
-0.912180E+02 
-0. 9xa5loE + 02 
-0.9 15J30E+02 
-0.91o711E+02 
-0. 9 1 8a95 E+ 02 
-0.9107 1 lE+02 
-0.92a094E+02 
-0.92a787E+02 
-0.928840E+02 

— 0. 9oial oE+Oi 
-0.96O298E+02 
-C. 9 + U87oE+02 
-0.946171E+02 
—0.9 5.1O34E+02 

— 0. 9 a 9a 4 7 E+Ox. 

— 0.968035E + 02 
-0.978J7aE + 0l 
-0.969995E+02 
-0. 100417E+06 
-0. 102iOaE+03 
— 0.104l01E+o3 
-O.10O4 50E+03 
-0. 109i8l£+0a 

— 0. 1 lioOOt+Ua 
-0.1157a4E+03 
— 0.119567t+0a 
-0.12854bE+06 
-0. xi7o00E+U6 
-u. 13 1594E+0a 
-0.185427E+03 
-0. i390aaE + 03 

— 0.1 424 O'+E+oa 
-0.i45aoiE+0a 

— 0. I4o567c+03 
-0.l51495E+0a 
-0. ia44i9E+06 

— 0. 157444E+0J 
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Table 1. (Contd.) 



0 .464157E + C1 
0 .3i.4i 74E + C1 
0.569579E+01 
0.63U955E+Cx 
0 .095 943 £+Ui 
0. 77426iE+0i 
0 .55769^£+0i. 
0.95Qli.6£+0i 
0 .105230E + 02 
0 .il559iE+02 
0-129 155E+02 
0.x43 07w:E + 02 
0 .l554d9E +02 
0.175567E + 02 
0 .x944d5E+02 
0.215443E+02 
0.238658E+C2 
0-2£>4375 E+02 
0.2 92 863E + 02 
0 .324421E + C2 
0 .3593OQE+02 
0.398i0o£+C2 
0 .^41004E+02 
0.488 526E+02 
0.541 167E+02 
0 .59;y483 E+02 
0 .664 OdOE + 02 
0 .7356 j9E+02 
0 .814910.E+02 
0.902 722E + 02 
0 .999996E + 02 



— 0. 48o013E + 01 
-0.636283E+01 
-0.768934E+01 
-0. 9449o6E + 0l 
-0. 110539E+02 
-0. 127116E + U2 
-0. 1443l3E+0^ 
-0. 162196E+02 

- 0. 180814E + 02 
-0.20O195E+02 
-0.2203^0E+02 
-0. 24x236E + 0^ 
-0.262850E + 02' 
-0. 285134E + 02 
-0.308063E+02 
-0.3314O6E+02 
-0.355431E+02 
-U.379807E+02 
-0. 4U4556E+02 
-0. 4^9o 26£ + Oz 
-0. 454962E+02 
-0. 480529E + 02 
-0. 50O288E + 02 
-0.532208E+02 
-0. 55826UE + 0^ 
-0. 5O4423E+02 
-0. 610677E + 02 
-0. 6o7005E + 02 
-0.66J395E+02 
-0. 689835E + 02 
-0. 716317E+02 



-0.1o0608E+06 

— 0.1639 74E + Oo 
— 0.io7576£+03 
— 0.i.7143^E + 0o 
-0. 1 7554ut + 0o 
-0.179o83£+03 

— 0. i84‘t j3E + 03 
-0.189136E+03 
-0. 19393o£+08 
-C.198782E+03 
— 0.20jo12c+U3 

— O.^OSoo9E + Uo 
-0.213006E+0O 
-O.rii 74O9E+03 
-0.221735E+06 
-0.225776E+03 
-0.2 z9574E+Oj 
- 0.2o3121E+03 
— 0.2io416t+03 
-0.^694 o^E+0j 

— 0.24<:^6 7E + 0o 
-0.244b^0£+03 
-0.247195E+0O 
— 0 .2490^4 E+ 03 
-0.^:5 x602E + 0o 
- 0.2530 83E+0^ 
-0 . 2 34700 1+03 
-0.^3olo8E+0o 
-0.237498E+03 
-0.258703E+03 
-0.259793E+O3 
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Table 2. Tabular Output of Example 4.2 



FRES 



U.xOQOOOE+00 
0 .ii0776E+OG 
0.L22 7 LZE + CO 
0 .J.JP935E+00 
0. 1503d4E+0U 
0.166810E+CG 
0.i.847d5E+00 
0.X04697E+G0 
G •34.0 7 34 E + OG 
0 .25ii88E+0G 
0.278356E+CG 
G *3 00 240 8 + 00 
0. J41454E+CG 
0.378248E+0G 
0 .41::^ GOO E+GG 
0.+64158E+CG 
0 •3141 75c + 0G 
0*369500.E+OG 
G.03G956E+CU 
0 •O90946E + GG 
0* 7742o3E+0G 
0*857 OS5E+00 
0 *9501178+00 
0*lG5250E+0l 
0 *lj.o591E + 0x 
0*1291558+01 
0*143 0 728 + 01 
0 *158409, E+Oi 
0*175567E+0l 
0*i.94485E + 01 
0 *215443 8+01 
0 * 238 65 8E + Cl 
0 *264375E + 0i 
O*292064£+Ol 
0*324421E+C1 
0 *359301 c+Oi 
0*o90l06.E + 01 
0 *441 004E + 0i. 
0 *4oo526E+0i 
0*5411688+01 
0 *599483E + 0i 
0 * 664 081 C + 01 
0*735640E+C1 
0 *ol4911E+0i 
0*9027238+01 
0.999999E+0X 
0 *llQ775E+02 
0*122712E+C2 
0 * J.J5935E + 02 
0*150583E + 02 
0*166809E+02 
0 *104 7o4E+02 
0*204696E+02 
0 *426754E + 04 
0*251 loOE+02 
0.27O255E+02 
0 *3082398+02 
0*3414548+02 
0 *o78247 8 + 03 
0 *4190078+02 



MAuiMl TUJE 



0* 590345w + 03 
0* 50V089E + 02 
0* 5797548 + 03 
0.5703258+02 
0*5607838+03 
0* 551 1 068 + 03 
0*5412708+02 
0* 5312498 + 02 
0* 52l0l3£ + 03 
0* 5105298 + 04 
u* 4997658+02 
0* 48boo7E + 02 
0*4772668+02 
0* 4o54748 + 03 
0* H532098 + 02 
0* 440O98E + 02 
0* 4376978 + 03 
0*4142928+03 
0* 40G5038 + 02 
0*30635oE+02 
0*3710918+03 
0* 3574558 + 03 
0*3422018+02 
0* ^270908 + 02 
0* 31J.081E + 03 
0*29oo408+03 
0* 2814328 + 03 
0* 2 oo 320E + 02 
0* 3513678+03 
0* 23oo3l8 + 02 
0*2221648+02 
0* 20d014E + 02 
0* 1942178+02 
0* 1008018+02 
0* 1677838+04 
G* 1551o68+03 
0* 1429478+03 
0* 1311078 + 03 
0* 1196248+03 
0 * 108+668 + 02 
0. 9759948+01 
0* 8698608+01 
0* 7658778 + 01 
0 * 66365o 8 + 01 
0* 56202OE + 01 
0* 4630068 + Ol 
0* 3630538 + 01 
0. 2650118 + 01 
0* 1661308 + 01 
0*6606568+00 
-0*361 j 96=+0U 
-0* 1343528+01 
-0* 2368478 + 01 
-0* 34l315E + 0i. 
-0*44605 oE+01 
~0. 557+838 + 01 
-0*6699928+01 
-0* 7859608 + 01 
-0* 9057388 + 01 
-0. 1039618+02 



PHASE 



- 0 * 995141 8+03 

— 0*1 0 05 03 8+ 03 
-0*1 Ol58e 8 + 03 
-0*1027698+03 
—0*10 40 5o 8+03 
-0* 1054578+03 
-0* 106966 8+06 
—0*1 085948+03 
-0*il033lE + 03 
-0*1121748+03 

— 0*ll4ilDE+03 
-0*1 161368+03 
-0*1164268+03 
-C* 1303568+03 
-0*1244918+03 
-0*1 2 46 06 8+03 
-0* 1466638+03 
-0*1266238+06 
-0*1304498+03 
-0* 132i 068+06 
— 0* 1365ouE+06 
-0*1647858+06 
-0*1357568+03 
—0*1 36+618+03 

— 0*l3 6084t+u3 
-0*1370228+03 

— 0*1668778+06 
-C* 1364578+03 
-0* 1 3 57 ?6 8+06 

— 0* 16+8538+03 
-O*l337l6t+06 
-0*1323958+03 
-0* 1309278+06 
-0*12965 iE+03 
-0.1377098+03 
-0* 1360438+03 
-0*1243938+06 
-0*1437978+06 
-0*1213898+03 
-0*1190998+06 
-0*1166528+03 
-0*1175698+06 

— 0*1 166 66 t+Oj 
-0*1159558+06 
-0* 1 1 5+45 8+06 
-0* 1161458 + 06 

— 0 * 1 1 50 60 8 + 06 
- 0 * 1 151928+06 

— 0*11 55 478+03 
-0* 1 161358+03 
- 0 * 1169378+06 
-0*1 179538 + 03 
- 0 * 1 194618+03 

— 0* X 4 Oo 66 8+ 03 
-0 *1333468+06 
-0* 1 244 19 8+03 
- 0 * 1262818+06 
-0*1 265108+03 
-0*1608838 + 06 
-0*1366768+06 
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Table 2. (Contd.) 



0.4o4137E+ C2 
0 .514174E+C2 
0.569579E+02 
0.630955E+C2 
0.69a 945 E+02 
0.774261E+02 
0 .857692E + 02 
0.950ii6£ + 02 
0.i05250E+CJ 
0 .iXo591E + 0i 
0.i29155E+03 
0.143C72E+CJ 
0 .158489E + C3 
0.A75567E+03 
0 .194485E + C3 
0 .215443 £+03 
0.238658E+C3 
0 .Z64375E + C3 
0.292 8 63 E+03 
0.324421E+C3 
0.359380E+C3 
0.3S6iC5E+Ca 
0.4410C4E+C3 
0 .4885^5E+03 
0 . 541 1 67E + 0 J 
0 .59948ZE + C3 
0 .664 080 E + 0 J 
0.735639E+C3 
0 .8149x0E+03 
0.902 72iE+C3 
0.9999S8E+63 



-0. ii5780E + 02 
-0. 129043E+02 
-0. i42758E+02 
-0. 166905E + 02 
-0. i7i484£+02 
— 0. ldt>46oE + 02 
~0. z0182^E+0^ 
-0.2i7522£+02 
-0. 233528E + 02 
-C..:498C6£ + 02 
~0. 266322E+02 
-C. 283043E + 02 
-0.29994l£+02 
-0.3io986E+02 
-0.334i61E+0^ 
-0.351439£+02 
“0. 368a05£ + 02 
-0. 38o245E + 0z 
-0. 403745£+02 
-0. 421296E+02 
-0. 438888£ + 02 
— O.456513£+0^ 
-0. 47416 6E + 02 
-0.491842£+02 
-0. 509536£+02 
-0.5Z7247E+02 
-0. 544969£+0^ 
-0. 56 z7C2£+02 
-C. 58044JE+02 
-0.596190£+02 
-C. 615944£+0^ 



~0.i35960£+03 
-0. i3o58i£+03 
-0.141232£+0a 
-0.i4387lt+0^ 

— 0. i 4646 7£ + 0a 
-0.148993£+03 
-0. l5i4^ 7£+0a 
-C.153749E+U3 
-C.155948E+03 
-C.x58014£+0ji 
-0.159944E+0J 
~0.1 o176o£+ 03 
-C.163392E+03 
-0.164917E+03 
~C. loojloE+Oj 
-C.167597E+03 

— 0.168765 E+0 3 
-C.l69830£+0a 
-0.170798E+03 
-C.171o78E+0a 
-C.172476E+03 
-0 .17 J199E+03 

— 0. X 7 j8 54 E+ Oj 
-0. 174447E+0J 

— 0 .1 7 4984 t+U 3 

— C.x7 5469E+0a 
-0.175907E+OJ 
-0.1 76a04E+03 

— 0.17 6o 6 2 £ + 0 3 

— 0.17 6966 £+ 0 J 
-C. 1 77z 78E+0 j 
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3. 



Example 4.3 : Velocity Feedback 



a. Problem Statement 

A simple second-order servo is to be 
compensated with tachometer feedback. The forward transfer 
function is 

100.0 

G(S) = 

S(S + 1) 

and the tachometer transfer function is K^S . The tachometer 
is fed back around all of the forward gain. Using Bode 
diagram methods, set to provide = 1.3. 

b. Solution 

Bode Plot for this system is shown on figure 
6. The system has a phase margin of about 6 degrees. For 
Mpw = 1.3, a phase margin of 45 degrees is required. A rough 
graphical design on the uncompensated Bode Plot gives 

1/H = 12.0/S or H = 0.08S 

The Bode Plot for the compensated system is on figure 7, 
showing a phase margin of 50 degrees. 

The close loop frequency response of this 
example, drawn using BODPLT is shown on figure 8. 
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Table 3. Tabular Output of Example 4.3 



FREw 



O.iOOOCOE-Ci 
0 * xiJ646E-Ci. 
0.129155E-Ci 
0.i467aOE-Cl 
0 .166 6 10 E~0 i. 
0.189573E-C1 
0 .2J.5443E-C1 
0 .244843 E-01 
0.27825oE-Ci 
0 .J162^7E-C1 
0.359o8iE-Cl 
0. 408 42:>E- Cl 
0 .464159E-C1 
0.527499E-C1 
0 .599483E-C1 
0 .631291E-01 
0.774263E-C1 
0 .879921E-Ca 
0. 999SS8E-C1 
0.113646E+CC 
0.129135E+00 
0 . 146 780E + CU 
0 • 166 8 ID E 4CG 
0 .189573 E+CO 
0.215443E+CC 
0 .^44843£ + CC 
0.278255E4CG 
0.31o227E+C0 
0 .359381 E+00 
0.408423E+CO 
0.464157E+CC 
0.527498E+CO 
0.599483E+CC 
0 .ot)i 291 E + CC 
0.7742c2E+CO 
0.879919E+CC 
0 .999996 E+00 
0 . 1 13 646E + C 1 
0 .129 i55E + Cu 
0.146779E+01 
0.1668C9E+C1 
0.i89573E+Cl 
0. 215-443 E + Cl 
0.244843E4C1 
0.276255E+C1 
0.316227E+C1 
0.O59381E+C1 
0. 408 4.^3 E •►Ox 
0.464X57E+C1 
0 .527498E + C1 
0.399482£-^Cl 
0.681290E+01 
0 .774261E + 01 
0.879919E-*-Ci 
0 .999996E + C1 
0.113646E+C2 
0.129155E+C3 
0.146779E+02 
0 .166809E + C2 
0.X89573E+C2 



MAGNITUDE PHASE 



0. 606303E + 02 
0. 595a9^E + 0x: 
0.584081E+0^ 
0. 5 72 9 70E + 0.:: 
0.561859E+02 
0. 550747E+02 
C. 539636E + C2 
0. 528525E ►■02 
0. 517414E-*-02 
0.5C63C3E + 0^ 
0.49519tE■^02 
0. 484080E-*- 02 
0. 472969E■^02 
0. 461858£-*-02 
0. 450746E + 02 
0.439635E+02 
0. 428523E + 02 
0.4174xlE-*-02 
0. 4 062 98 E ■►02 
0.395186E^^0.d 
0.384073E^^02 
0. 37z959E^^02 
0.36l845E-^02 
0. 350730E-^02 
0. 339613E^^02 
0. 328495E^^02 
0.3l7j7t>E + 0Z 
O. 306253E^^0.^ 
U,295127E+02 
0.283997E + 02 
0. ^72862E^^02 
0.261720E-^0z 
0 . 2503t8E + 02 
0,2^9404E+02 
U.228226E+U2 
0. ^a 7 J28E + Oz 
0. 205604E+02 
0. 194349E^^02 
0. 183252E^^0.^ 
0.171902E+0Z 
0. 160484E^^G2 
0. 14i3980E-^02 
0.137367E^^02 
0. lzo615E+ 02 
0.1 13t.91E-^02 
0. i01552E^^Oz 
0. 891480£^^01 
0 . 764.C22E401 
0. 63j 108E-^01 
0.497453E-^0l 
0.35o567E-^0i 
C. 2C9 793E4 0.t 
0.5657 C6E-^0u 
-0. 1035C6E-^01 
-0..<:70o49E-^0i 
-0.444846E-^01 
-0. 6 .c: 5853E-^C1 
— 0.6 132 4 OE-^ Ox 
-0. 100643E^^02 
-C. lz0476E + 0z 



-C.90O6 04E-^02 
-C,90Q6aoE+02 
-C,9007ii^E+0z 
-0.90Oo9zE-^0x: 
-C.901016E^^02 
-0.90ll56E^^02 
-C.901j1oE^^02 
-C.901497E+U2 
-C.901703E^^02 
-C.901936E^^02 
-ij.902202E+02 

— 0. 9 0 25 04 E-^Oz 
-C.902848£-^02 
-0.90J238£■^02 
-0.9036 8 zE-^Oz 
-C.904185E-^02 
-C.90475oE^^ 02 
-C.905409£-^02 
-0.90ol49E^^02 
-0.90o989E-^0z 
-C.907945E-^02 
-0.909u30E+02 

— C . 9 X 02 o 4 E •► 0 z 
— C.9il666E+0Z 
-C.913Z59E^^02 
-C.915069E-^0 z 
- 0.9l71ZoE^^0z 
-0.91 9466 E-^Oz 
-C.922118E-^0z 
-0.925134E-^02 
-C.9.t85o0E-^0Z 
-0.93Z451E^^U2 
-C.936870E+0Z 
-C.94188o£+02 
-0.947579E-^02 
-C.954037E^^0Z 

— C.96 1360E+U2 
-0.9o96 5 7E-^02 
-C.97905zE-^ 0Z 
-0.96 9o76t-^02 
-C. 100167E+03 
-C. 10x520E-^03 
-0. 106042E-^03 
- C. xO 47 48 E-^06 

— C. X 0666ot + 03 
-O.x 087 78 E-^03 
-C.1 x1iz7E-^06 
-0 . 1 X6708E+06 

— 0.11 6522 E "^06 
-C.119560E+06 
-0.1z2804E-^03 

— C. x2ozZ4E-^U6 
-C.129777E^^03 

— O . 1 334-1 3 E"^03 
-C.137075£^^03 
-0.140704E^^06 
-O. i44z4zE-^03 

— C.147639E"^03 
-0.150867E-^03 
-C.x5686oE-^03 
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Table 3. 



(Contd . ) 



0.2 15443 E 
0.Z44843E 
0 .278255E 
0 .316226 E 
0.359380E 
0 .4U8423E 
0.4o4157E 
0.527497E 
0. 599482 E 
0.681290E 
0 .774261E 
0 .879919E 
0.999S95E 
0 .li.3c46E 
0.129154E 
0.146779E 
0 .166 809 E 
0 .i.89573E 
0 .215^4jE 
0 .244 843 E 
0.278255E 
0 .Jlo^26E 
0.J59379E 
0.^08422E 
0.464157 E 
0.5274S7E 
0-599481E 
0 .681290E 
0.774260E 
0 .8799a8E 
0.999997E 



*02 


-0.140753E+0^ 


~C.i5 6650E+0j 


+ C2 


-0. 161408E+0^ 


-0.159199E+03 


+ C^ 


~ 0 . i. 8z3 7 6E + 02 


-0. 1 6153 7E + 0j 


+ 02 


-0. 203o02E+02 


-C.io3609E+03 


+ C2 


-0. 225037E + 02 


-0.1o5489E+03 


+0^ 


—0. 24o639E+C2 


-0.1o7170E+03 


+ C2 


-0.26O375E+02 


“0.1006 68 E +03 


+ C2 


-0. 290217E+0^ 


“0. io9999E+03 


*02 


-0.312143E+02 


-0. 171179E+03 


+ C2 


~0. 3341j4£+02 


-0.172224E+03 


+ 0z 


-C.35oi77E+02 


-0.173i48E+C3 


+02 


-0.376259E+02 


-0.1739648+03 


+ C2 


-0. 400373E+02 


— 0.17 4o84 E+03 


+ CJ 


-0. 4^^5llE+02 


-C.175319E+03 


*C3 


-0.444O68E+02 


-0.175o79E+0o 


+ C3 


— 0. 466840E + 02 


“0.17 03 7oE+U3 


+ 03 


-0.489023E+02 


“0.17 68 06 £+ 03 


+ C3 


-0. 511215E+02 


“0.177189E+03 


+ C3 


-0. 533413E+0.: 


“0.l7753oE+03 


+ 03 


-0.555617E + 0.i 


“C.177O22E+03 


+ C3 


-0. 577825E + 02 


“0.17S083E+03 


+ C3 


-0.600037E+02 


“0. 17O313E+03 


+ C3 


-0.622250E+02 


“0.178515E+03 


+ C3 


-0. 644466E+0Z 


“0. A 78o9oE+0o 


+ 03 


— 0.66o683E+02 


“0. 178O50E+03 


+ 03 


-0. 688902E+02 


“0.178983E+03 


+ 03 


-0. 7l3l^lE + 02 


“C.x79i09E + 03 


+03 


-0. 73334 lE + 02 


“C.l7 9215E + 03 


+ 03 


— 0. 75556i.E + 02 


“C.179309E+0O 


+ 03 


- 0 . iiiiazB+oz 


“0. 17939iE+0o 


+ 03 


-0.8U0003E+02 


“0 .1 7 946 5 E+03 
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Figure 8. Closed Loop Response of Example 



V, 



CONCLUSION AND RECOMJ-IENDATIONS 



A. CONCLUSION 

The objective of this thesis was to develop an inter- 
active user oriented computer program which would aid in 
solving control engineering problems using the Bode method 
of design. The presented program proves that frequency 
domain design of control systems using the digital computer 
as an aid is not only feasible but highly desirable. The 
results obtained are readily interpretable and provide good 
and meaningful insight into the problem. 

The results obtained during the investigation of the 
program performance show that a complicated but well formu- 
lated problem can be solved with ease, and the solution is 
obtained with speed, accuracy and precision. 

The entire program is less than 1000 lines, with a total 
of 9 subroutines. Every effort has been made to keep the 
program simple yet unambigous, so that the user has to 
invest very little time learning how to use it. Effort has 
also been made to minimize the use of the computer CPU time. 
However, expenditure of CPU time is to a large measure 
dependent upon: 

(1) the order of the system. 

(2) the number of iterations used in reaching a 
solution. 
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(3) the type and order of the compensator (s ) used. 



B. RECOMMENDATIONS 

The program as presented in this thesis seems to be able 
to adequately satisfy the usual needs in a control system 
design problem. A number of useful extensions to the work 
developed in this thesis can be carried out. These are 
briefly discussed. 

1. Curve Fitting 

Although not specifically worked on in this thesis, 
the program can be used quite effectively for curve fitting 
purposes. This was demonstrated in the initial stages of the 
development of this program. The procedure is, by its very 
nature, iterative and therefore time consuming and cumbersome. 
However, the entire algorithm can be automated using a 
minimization subroutine, with the program outputting a 
polynomial to fit a given curve over a specified range of 
the independent variable. 

2 . Computer Selection of Compensators 

A suitable minimization routine such as Box PLX 
can be incorporated into the program which could select the 
best possible location of poles and zeroes to meet given 
performance specifications. This would automate the entire 
Bode design procedure, the user then having the option of 
only specifying the type of compensation, i.e., cascade or 
feedback. It may however be pointed out that minimization 
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routines by the very nature of their operation are very time 
consuming and wasteful of CPU resources. 

3 . Root Locus 

Most of the subroutines developed in this thesis 
are very general in nature and can be adapted quite easily to 
develop a similar interactive program for Root Locus plots. 

4 . Integrated Control System Design 

No meaningful design of control systems is complete 
without finally analyzing its time domain performance. It is 
therefore considered highly desirable to incorporate into 
this program, an interactive Root Locus design procedure and 
then a time domain analysis of the compensated system. The 
entire package would then be an excellent teaching aid for 
control system design. 
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APPENDIX 
PROGRAM LISTING 



OOOCJOOOOOOOOOOOOOO 0 00000^0 0:30C3 0000 03 DOOOOOOOOOOO 

-n or^oo oo>o hcnj o^n 04- r^Oh- X) 
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